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DISTANCE PROTECTIVE RELAY USING 
A PROGRAMMABLE THERMAL MODEL FOR THERMAL PROTECTION 
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Technical Field 

This invention relates generally to power line 
thermal protection, and more specifically concerns such thermal 
protection implemented with a first order thermal model 
10 emulating the temperature of the power line. 

Background of the Invention 

Thermal protection for power lines is accomplished by 
monitoring the temperature of the wire conductor and generating 

15 a trip signal which opens a circuit breaker when the conductor 
temperature becomes greater than the maximum allowed temperature 
for the conductor. The temperature of the conductor is not 
measured directly, but is obtained by using a first order 
thermal model for the conductor. This involves using a heat 

20 balance equation (heat input minus heat losses) of a 1000 foot 
section of conductor. 

The 1000 foot section of the line is assumed to have 
the highest temperature and to be exposed to the maximum solar 
radiation. The heat input portion in the equation is mainly due 

25 to heat dissipated in the conductor resistance and the solar 
heat gain. The heat loss portion is due primarily to convection 
and radiation. In general, the equations for establishing such 
a thermal model are well known, and the procedures for obtaining 
the temperatures of overhead conductors using a thermal model 

30 are set forth in IEEE Standard 738-1193, titled "IEEE Standard 
for Calculating the Current -Temperature Relationship of Bare 
Overhead Conductors" . 

A power line thermal model has been used previously 
in protective relays, although the thermal model equations are 

3 5 preestablished and the parameters thereof are determined by the 
manufacturer. The user has no control over the thermal model or 
its operation. It is an automatic process. This can lead to 
inaccurate results in some cases. Accordingly, it is desirable 
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that the end user of a protective relay have some capability to 
itself establish the thermal model parameters for a more 
accurate and fast determination of power line temperature. 

Summary of the Invention 

Accordingly, the present invention is a protective 
relay system for power line thermal protection, using 
programmable logic present within the relay, comprising: a 
protective relay for power lines, which includes a programmable 
logic capability by which the end user of the protective relay 
can enter operational settings which are then used by the relay 
in carrying out its thermal protection functions; a set of 
stored thermal model equations which- emulate the temperature of 
a power line conductor, based on a plurality of individual 
setting values which are enterable into the relay by the end 
user, and wherein the logic implements the entered setting 
values into the thermal model equations which produce an 
emulated temperature of the conductor; and means for providing 
an indication when the temperature of the conductor exceeds a 
preselected value. 

Brief Description of the Drawings 

Figure 1 is a block diagram showing a general 
implementation of a thermal model for determining temperature of 
25 a power line conductor. 

Figure 2 is a block diagram of a particular example 
of a thermal model implementation. 

Figure 3 is a block diagram of the processing used by 
a thermal model to produce a temperature value. 
3 0 Figure 4 is an example of a program code for a 

thermal model . 
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Best Mode for Carrying Out the Invention 

As indicated above, use is made of first order 
3 5 thermal models to determine the temperature of electrical 
conductors, such as for instance the temperature in various 
parts of an induction motor. These principles are set out in a 
paper by S. A. Zocholl and Gabriel Benmouyal entitled "Using 
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Thermal Limit Curves to Define Thermal Models for Induction 
Motors" , Proceedings of the 28th Annual Western Protective Relay 
Conference, Spokane, Washington, October 2001. The first order 
thermal model for protection of a power line is expressed as 
5 first order differential equation: 

d(TC) 

Power Supplied (P) -Losses = (1) 

dt 

The equations below use the following variables: 

P s= heat power input to the conductor (kW/kft) 

THC = conductor heat capacity (kJ/°Ckft) 
10 TRA = thermal resistance to ambient (°Ckft/kW) 

TC = estimated conductor temperature (°C) 

TA = ambient temperature (°C) 

TI = conductor initial temperature (°C) 

I = conductor current (A RMS) 

15 rac = AC conductor resistance at 25°C (Q/kft) 

rdelt = temperature coefficient of AC resistance (Q /°Ckf t) 

Qsun = heat power input from the sun (kW/kft) 

Qradiated = radiation heat losses (kW/kft) 

Qconvected = convection heat losses (kW/kft) 
20 At = processing interval or integration step (s) 

Basically, the heating power supplied to the conductor can be 
expressed as: 

P = I 2 • (rac) + (TC - 25) • rdelt) + Q sun (2) 
25 Further, the heat losses from the conductor are expressed in the 
following equation: 

TC-TA 0 

Jiftj^ Eradiated Qconvected v ' 

The conductor temperature differential equation is then: 

THC = I 2 . {rac + (TC - 25) . rdelt) + Q sm - ( 4 ) 

at 1KA 

30 From Equation (4) , the solution for the conductor temperature is 

obtained as follows: 

TC(t) = if ^^acHTC-25).rdeltWQ sun _ TC-TA V 
w l\ THC TRA •THC) 
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The equation is solvable numerically by computing the 
conductor temperature in successive time increments ATC, to 
determine the conductor temperature . 



ATC = 



^ I 2 • (rac + (TC - 25) • rdelt) + Q sun TC-TA ^ 



AT (6) 



THC TRA • THC ) 

5 The conductor temperature is determined successively 

as follows: 

TC new = ATC + TC old (7) 

The thermal model is illustrated in Figure 1, where TTH is the 
maximum allowable temperature threshold, above which the circuit 
10 breaker for the power line can be tripped, while TTL is an alarm 
threshold. When Q sun is assumed to be a constant, and when t is 
infinite, the conductor temperature is represented by the 
following equation: 



= / T 
\ (rac + 



TTH -TA- Q sun • TRA , % 

^Csun (8) 



(TTH - 25) • rdelt) • TRA 
15 Referring to Figure 1, -the heat power into the 

conductor is represented by source 12, while THC (conductor heat 
thermal capacity) is represented by capacitor 14, TRA (thermal 
resistance to the ambient temperature value) is represented by 
resistance 16 and TA (ambient temperature) is represented by 
20 battery 18. The result of the thermal model is an estimated 
conductor temperature TC at output point 20. The temperature 
value at output point 20 is compared to the value of TTL (input 
21) in comparator 22 to produce an alarm, and compared with the 
value TTH (input 23) in comparator 24 to produce a trip signal. 
25 Several specific values in the above equations must 

be determined in order to solve the equations and complete the 
first order differential equation ... The value of Q sun is 
determined by the following equation: 

Qsun = SAC • DIA • SIR (9) 
3 0 The variables in the determination of Q sun are set forth below. 

d = solar declination (degrees) 

n = days of the year (days) 

Hrs = local time in hours (hours 

LocT = local time (degrees) 
35 Ws = solar time (degrees) 
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solar zenith angle (degrees) 


LSTD 




standard time meridian (degrees west longitude) 


LON 




conductor longitude (degrees west longitude) 


LAT 




conductor latitude (degrees north latitude) 


SIR 




solar incident radiation (kW/in) 


Qsun 




solar heating (kW/kft) 



10 Referring to Equation 9, the value of SIR is 

typically provided by a look-up table as a function of the solar 
zenith angle. The solar zenith angle is determined indirectly 
by its cosine as a function of the solar declination d, the 
conductor latitude LAT and the solar time Ws as in the following 

15 equation: 

cos(z) = sin(rf) • sin(LA T) + cos(c/) • cos(LA T) • cos(Ws) (10) 
The value of the solar incident radiation (SIR) can then be 
computed using the following equation: 

SIR = a 5 • cos(z) 5 + a 4 • cos(z) 4 + a 3 • cos(z) 3 + a 2 • cos(z) 2 + a x • cos(z) 1 + a 0 ( 11 ) 
20 Power n of cos(z) can be determined by multiplying the number 
n times by itself. The solar time Ws is provided as a function 
of the local time LocT, the longitude value LON and the standard 
time meridian LSTD as follows: 

Ws = LocT + (LON - LSTD) ( 12 ) 

25 LocT is determined by a function of the local time expressed in 
hours, as set forth in the following equation: 

LocT = (Hrs - 12) • (-15) ( 13 ) 

The solar declination d can be expressed as a function of the 
day of the year n, as follows: 

30 d = 23 .45 • sin( 284 + " • 360) (14) 

365 

The local time and hours in the day of the year n would be both 
available within the relay as THR and DDOY. The conductor 
latitude LAT and longitude LON, along with the standard time 
meridian LSTD, must be entered as set-tings. 
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Another determination which must be made is the 
conductor thermal resistance TRA, which is determined by 
Equation 15 below: 

TRA = TC-TA (is) 

Qradiation Qconvected 

5 and is equal to the temperature difference between the conductor 
and the surrounding air, so that the heat transferred to the 
conductor is one unit of power. 

The radiated heat losses can be computed as follows: 

Q rad = S • E • A*(KC 4 - KA 4 ) (16) 
10 with the following variables: 

S = Stephan-Boltzman constant (W/°k 4 ft 2 ) 
= 0.527e-8 

E = thermal emissivity constant (unitless) 
= 0.23 for new conductor 
15 =0.91 for blackend conductor 

A = area of circumscribing cylinder (ft 2 ) 
KC = conductor temperature (°k) 
KA = ambient temperature (°k) 
After adjusting the units, 



25 



where DIA (the conductor diameter) is in inches and the units of 
Q rad are in kW/kft . 

The value of the convected heat losses (Qconvected) is 
determined by the following equation: 

,0.52 



Qcon* = [1.01 + 0.371 ^^*^ . }.K. (TC-TA) (18) 



where the variables are: 

R = air density (lb/ft 3 ) 
V = air velocity (ft/hr) 
H = absolute viscosity (Ib/hr ft) 
30 K thermal conductivity (W/ft 2 °C) 



The value of thermal capacity (THC) is determined by 
multiplying the number of pounds of aluminum (WA) by the 



7 

aluminum specific-heat and adding the number of pounds of steel 
(WS) multiplied by the steel specific heat, as in the following 
equation: 

EM* 428.8 + ETC* 204.9 (19) 
100 

5 The above equations, as previously indicated, used to 

determine a thermal model, are known. However, previously a 
manufacturer of a protective relay, using such a thermal model 
for temperature determination of conductors, would directly 
provide all the specific variables at the factory and completely 

10 set up the thermal model in the relay for the customer. The 
customer thus has no ability to affect the thermal model in this 
arrangement. Frequently, however, this would be less than 
optimum, since the customer's use/location would likely increase 
the accuracy of the thermal model . 

15 in the present invention, a protective relay includes 

a logic capability by which the user can directly program the 
thermal model equations in the relay. In order to provide the 
information necessary for determination of the thermal model, 
the following settings are made by the user. The first group of 

20 settings is the solar model settings. 

Solar Model Settings 

DSH = Default solar heating (kW/kft) 
SAC = Solar absorption coefficient (unitless) 
25 DIA = Conductor diameter (in) 

LSTD = Longitude of time standard (degrees west longitude) 
LON = Longitude of conductor (degrees west longitude) 
LAT = Latitude of conductor (degrees north latitude) 

The diameter of the conductor is corrected, both according to 
3 0 elevation and atmosphere, in accordance with the following 

equation: 

DIA cor =k x *k 2 *DIA (20) 
Values of k x and k 2 are shown in Tables I and II provided below. 
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Table I Table II 

Elevation (ft) kl Atmosphere k2 

0 1.00 Clean 1.00 

5000 1.15 Industrial 0.82 

10000 1.25 
15000 1.30 

The longitude of the time standard is entered 



according to Table III 



Table III 



10 Time Zone Meridian 

Eastern 75°W 
Central 90°W 
Mountain 105°W 
Pacific 120°W 

15 Next, the thermal model settings as set forth below 

must be entered. 

rac = AC resistance at 25°C (Q/k) 

rdelt = temperature coefficient of AC resistance (Q/°C kf t) 
THC = thermal heat capacity (kJ/°C kft) 

2 0 TRA = thermal resistance to ambient (°C kft/kW) 

Next, the temperature settings must be entered. The 
estimated offset temperature EOT is "the temperature difference 
between the ambient temperature available and the estimated 
hottest ambient temperature along the line. It must be entered 
25 as a number other than zero if it is estimated that a difference 
in fact exists. If not, it is entered at zero. 
EAT = estimated ambient temperature (°C) 
EOT = estimated offset temperature (°C) 
TTH = high temperature threshold (°C) 

3 0 TTL = low temperature threshold (°C) 

TI = conductor initial temperature (°C) 

Lastly, the following logic switches must be set. 
THE = thermal trip enable (I/O) 
TSE = thermal sensor enable (I/O) 
35 SGE = solar generator enable (I/O) 
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THE has to be set to 1 if the line is to be tripped due to 
excessive temperature; TSE is set to 1 if an external sensor is 
available for the ambient temperature. SGE is set to 1 if Q SU n 
is computed from Equation (9) . It is set to zero if Q sun is 
introduced as equal to DSH. 

One example of a particular implementation is set 
forth below. A line conductor is made of Drake cable type with 
the following basic characteristics that can be extracted from 
tables : 



Diameter DIA=1.108 in (kl=kl=1.0) 

Thermal emissivity constant E=0.5 (unitless) 

AC resistance per mile at 25°C RL=0.117 Q/mile 

AC resistance per mile at 75°C RH=0.139 Q/mile 

Weight of aluminum 1000ft WA=750 Ib/kft 
of conductor 

Weight of steel of 1000ft WS-344 lb/kft 
of conductor 



The value of rac, the conductor resistance per 1000 ft, is 
determined from the value of RL: 

rac = — = = 0.022159Q/*/? (21) 

5.28 5.28 

The value of rdelt, the conductor temperature coefficient, is 
computed from: 

rdelt = RH -RL — = 0.139-0.117 = 0 0000833Q / o C ^ {22) 
(75 -25) #5.28 50*5.28 
The thermal heat capacity is determined by multiplying the 
number of pounds of aluminum WA by -the aluminum specific heat 
and adding the number of pounds of steel multiplied by the steel 
specific heat: 

= WA. 428.8 + MS .204.9 = 392 . 0S6kJ rC kft (23 ) 

1000 

TRA is computed using equations 15 through 18 and using a 
conductor temperature of 90°C and an ambient temperature of 
40°C. 

First, we compute the radiated power. From Equation 
17, we have: 
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30 



Q rod =0.l3&mE.D£A 



ATC , 4 _ , KA . A 
MOO 100' J 



(T^y-ir^y I < 24) 



or 



e rod = 0. 1 38 . 0.5 • 1 . 108 . - (^) 4 ^ = 5.937 kW/kft (25) 

The convected power is computed from Equation 18 with the 
5 following assumed air constants: 

R = Air density=0.0752 lb/ft 3 (20°C / sea level) 
V = Air velocity=2 ft/sec=7200 ft/hr 
H = Absolute viscosity-0 . 0439 lb/hr ft (20°C) 

K = Thermal conduct ivi ty=0 . 00784 watts/ft 2 °C(20°C) 
10 < 26 > 

Qcon* = [1 01 + a37 ( 1 ' 108 ^ Q439 * 7200 ) ]• 0.00784 #(90- 40) = 20.964W/^ 

Finally, we have for TRA: 

TRA= TC ~ TA = 90-40 = i m S59°Ckft/kW (27) 

Qra^Qcon, 5.937 + 20.964 

The time constant of the first order thermal model is provided 

15 then as: 

Thermal time constant = 
THOTRA=392 . 086»1 . 859=728. 3s=12 . 15 mn (minutes) (28) 

The practical meaning of this time constant is that if a step 
20 current is applied to the conductor, it will take 12.15 minutes 

for the conductor to reach 63% of the new steady-state 

temperature value . 

Assume that the line is, located in the Montreal 

(Canada) area. Also, assume the routine is started with the 
25 line open in summer time. The final set of settings that has to 

be introduced into the routine as constants as follows: 



DSH 


default solar heating 


= 0 


SAC 


solar absorption coefficient 


= 0.5 


DIA 


conductor diameter 


=1.108 


LSTD 


longitude of time standard 


=75.00 


LON 


longitude of conductor 


=73 .43 


LAT 


latitude of conductor 


=45.56 
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rac AC resistance at 25°C =0.022159 

rdelt temp, coeff. Of AC resistance =0.0000833 

THC thermal heat capacity =392.086 

TRA thermal resistance to ambient =1.859 

5 

EAT estimated ambient temperature =25 

EOT estimated offset temperature =0 

TH high temperature threshold =90 

TL low temperature threshold =80 

10 TI conductor initial temperature =25 

THE thermal trip enable =0 

TSE thermal sensor enable =1 

SGE solar generator enable =1 

15 

Figure 2 shows a block diagram for the above example 
to produce an alarm or a trip signal. The block diagram 
includes a thermal trip enable signal 25 to an AND gate 26 (with 
the output from the trip comparator) . The output of AND gate 2 6 

20 is applied to a timer 27, while the output of the alarm 
comparator 22 is applied to a timer 29. The timers 27 and 29 
are conditioning timers with pick-up and drop-out times of 10 
power signal cycles. 

As indicated above, settings entered by the user, 

25 such as in the example above, are implemented into the thermal 
model by a logic program which is present in the relay. The 
requirements of one such logic program as an example are set 
forth below. The program flow chart for programming the relay 
to implement the user-defined thermal model in the relay is 

3 0 shown in Figure 3. 

At the start of the program sequence, it is first 
determined whether or not it is the first processing sequence 
for the thermal model, as shown at block 30. If yes, the 
initial conductor temperature is provided, at block 32. If not, 

35 then if SGE (solar generator enable) is equal to one, as 
determined at block 34, then Q 8un is computed at block 36, while 
if not, it is set to the DSH (default) value, at block 38. The 
largest of the three phase currents (A, B, C phase) is then 
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chosen as the conductor current at block 40. The temperature is 
then calculated using the thermal model, at blocks 42 and 44. 
If the temperature is larger than the THL value (trip) , as 
determined at block 45, then the circuit breaker on the line is 
5 tripped. If the temperature value is not equal to or above the 
trip threshold value but above the alarm value, as determined at 
block 48, then an alarm bit is set at block 50. The sequence is 
then ended. The sequence repeats itself at specific intervals. 

As indicated above, applicants' invention includes 

10 the use of logic circuits and logic equations in the relay to 
accomplish a user programmable thermal model. In addition to 
the functional flow chart discussed above, one example of the 
logic coding to provide the programmable capability is shown in 
Figure 4. In this example, lines 1-19 introduce the power line 

15 thermal element settings, lines 20-39 provide the computation of 
Qsun# lines 40-42 provide the conductor initial temperature 
value, lines 43-46 determine the maximum phase current 
determination, lines 47-48 establish the ambient temperature, 
lines 49-52 establish the conductor temperature, lines 53-55 

2 0 determine the alarm output and lines 56-58 determine the trip 
output . 

The combination of the relay having the structure and 
capability of permitting the user to enter appropriate settings 
for a thermal model, along with the logic capability to 

25 implement those settings into the thermal model and then 
accomplish the calculation in accordance with the stored thermal 
model equations, results in a user-programmable thermal model 
system for protective relays. 

Although a preferred embodiment of the invention has 

30 been described for purposes of illustration, it should be 
understood that various changes, modifications and substitutions 
may be incorporated in the embodiment without departing from the 
spirit of the invention which is defined in the claims which 
follow. 

35 What is claimed is: 



